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1. INTRODUCTION 


Upon a surface of positive curvature there is a unique conjugate system for 
which the angle between the directions at any point is the minimum angle 
between conjugate directions at that point. ‘This system of lines is called 
characteristic lines or mean-conjugate net (Eisenhart, 1909, 1947). The object 
of the present paper is to obtain the equation of characteristic lines of a 
hypersurface in a Riemannian V,,,, and to discuss some properties of these 
curves. 


2. HYPERSURFACE IN Vn; 


Consider a hypersurface V,, of co-ordinates x‘, i = 1, 2,....n and metric 
ds* = gijdx'dx! ; (2.1) 
imbedded in a Vy,, of co-ordinates y*, a = 1, 2,....1 + 1, and metric 
ds? = a,gdy*dy?. (2.2) 
We then have the relation 
Sig = Gag Y*,iY*, j (2.3) 
where (,) followed by an index indicates covariant derivative with respect 


to the x with that index. If N*, a= 1, 2,....2+ 1, be the contravariant 
components of the unit normal to Vn, 


dag N®NF = 1, (2.4) 
¥* 3; ig = QN* (2.5) 
where (;) followed by an index indicates tensor derivative with respect to 


the x with that index and 2;; are the components of a symmetric covariant 
tensor of the second order symmetric in the indices i and j. 


3. CHARACTERISTIC LINES IN Vy 


Let ep, p = 1, 2,....n be the unit tangents to the lines of curvature 


inV,. If t, t be any two directions in Vn, we may express t, t as 
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t= ~ Cp, COS Op, (3.1) 
t= E ep, c0s ‘ (3.2) 
where @,, 0) are the angles which t and t make with the direction Cnr. 


The angle ¢ between t and t is given by 
cos $ = gij E ep,’ cos Oy F eg,) cos Og 
Pp qa 


= 2 cos 9, cos ‘.. (3.3) 
? 
since the congruences ep, (p = 1, 2,....m) are mutually orthogonal. 


Defining characteristic directions as a pair of conjugate (but non- 
asymptotic) directions, for which ¢ is extremum (so that cos ¢ must also be 


extremum), the characteristic directions are given by those values of 6p, Bp 
for which cos ¢ is extremum, subject to the relations 


2 cos? Oy = |, (3.4) 
= cos? 6, = 1, (3.5) 
and 
Qi; Z ep,’ cos Op J eg,4 cos Og = 
or 
Z Qij ep,iep,) cos Oy cos 8) = 0 


or 

= Z kp cos Op cos % = 0 (3-6) 
where k, is the principal curvature in the direction ep,. 

Applying the Lagrange’s method of undetermined multipliers, we have 


Z (sin % cos Opd0) + cos Oy sin F,d0y) 
? 
+ AZ sin Op cos Opd%p + XX sin G, cos Fpd0y 
? 
+E (kp sin 9) CoS Opd0p+Kp cos Oy sin 9,d8,,) = 0 (3-7) 


where A, A and » are the Lagrange’s multipliers. 











: 
: 
' 
i 
| 


the following set of equations for the determination of A, A, » and the @’s: 
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Equating to zero the coefficients of the differentials d0y, d0, we have 


sin 9p (cos Op = a cos 4p + pkp cos Oy) = 0 (3.8) 
sin 4p (cos Oy + Acos I) + pky cos Op) = 0 
p= 1, 2. ...™ 
But sin 6) # 0 for any value of p. For, if sin 0, = 0 for a fixed p, we have 
cos? 0, = 0,9 #p,q=1,2,....n in which case the direction ¢ must be 
a principal direction. Similarly sin 6, #0 for any value of p. 


Therefore the equations (3-8) reduce to 


cos 8, + A cos 9p + pk cos Op = Q, (3.9) 
Cos Ip + Acos Oy + pkyp cos) =0 |, 
pw BZ... 08 
Multiplying the first of the equations (3-9) by cos 4% and the second by 
cos A and summing up on p we have, 
cos¢+A=0 
cos¢ +A=0 \, 
by virtue of the equations (3 .3)—{3.-6). 


(3-10) 


From the equations (3.10) we have 
A=) (3-11) 
Using (3.11), we may write (3.9) as 
(1 + pkp) cos Op + Acos Oy = 
(1 + pkp) cos % + Acos 8 =0 (3.12) 


For a fixed p, say p=r, the solutions of (3.12) are either cos 6, = 0, 
cos 6, = 0 or else 


(1 + pk,)*? = A? (3-13) 
Similarly if for p = s, cos % #0, cos 9, #0, 
"(1+ ks)? = (3-14) 





















56 M. K. SINGAL AND RAM BEHARI 


There must exist at least two values of p, say r and s, such that cos Op, COS Op 
are different from zero; for otherwise, the directions t and t reduce to 
principal directions. 


Two different cases arise according as the principal curvatures are all 
unequal or there are sets of equal values of the principal curvatures. We 
shall discuss the two cases separately. 


(i) When all the curvatures are unequal.—From (3-13) and (3-14) it is 


found that there cannot be more than two values of p for which cos 4p, cos 8, 
may be different from zero. For, we cannot have more than two consistent 
equations of the form (3.13) for determining A and p. 


Thus, the characteristic directions are given by 


Cos Oy = 0, cos 9 = 0, p = 1, 2,....2, pFr,s (3.15) 
1 + pk, =A, 
weg \ (3.16) 
1+pks=—A 
since ky # ks. 
Using (3.16) we have from (3-12) 
cos 9, + cos 6, = 0 (3.17) 
cos 0, — cos 6; = 0 


Also using (3-15) we have from (3-4), (3-5) and (3.6) 
cos? 6, + cos? 6, = 1 
cos? 6, + cos? 0, = 1 (3.18) 
k, cos? 0, — ks cos? 6, = 0 


From the equations (3.18) we have 


tan? 6, = tan? J; = keke 3-19) 


tan? 6; = tan? 05 = ks/kr 
(3.15) and (3.19) define a pair of real directions provided k,, ks, be both of 


the same sign, showing that corresponding to every pair of principal curva- 
tures having the same sign we have a real pair of characteristic lines. 
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The angle between the characteristic lines as defined above is given by 
cos $ = cos 4p cos 4, 
= cos? 0, — cos? 6, 


= 2 cos? @, — 1 








= COs 20; = ke + ker (3 .20) 
The normal curvature in the characteristic direction t is 
k = k, cos? 6, + ks cos? 0, 
— 2krks 
af eS (3-21) 
Similarly the normal curvature in the characteristic direction t is 
k = ky cos? 0, + ks cos? 9, 
_ leks 


Thus the normal curvatures in the characteristic directions lying in the pencil 
determined by the principal directions e,, and es, are equal, each being equal 
to the harmonic mean between the principal curvatures k, and ks. 
(ii) When there are sets of equal values of the curvature at the point. 
Let kj =k,=.... =kyp Keg = Kegyy = .... = Kgiq-i, et. 
Following the same reasoning as in case (i) we find that corresponding to 
two distinct values of the principal curvatures, say k, and ks, there exist 
values of p such that 
(1 + pkp)? = »? (3.23) 
| cos 4) = 0, cos 9 = 0, (3.24) 
p=r-+1,r+2,....s—l,st+tqs+q+l,....n. 
In this case, (3.4), (3-5) and (3.6) give 
(cos? 6,++cos? @,.+....+ cos? 0,) 


) + (cos? 05 +....+ COS? 9519-3) = 

(cos? 6,-++cos? 6,+-....-+ cos? 6,) 
f + (cos? 6, +....+ cos? Gaal =1 tf = (3.25) 
- k, (cos 8, cos 6, +....-+ cos 9, cos 6,) 


+ ks (cos 6, cos 45 +.... 
— Cos 9si.q-1 cos Sea) = 0. 
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From (3.23) we have 


1 + pk, = A, ; 
+ iil eam , (3.26) 
giving 
2 k,—k 
tt peg Dee on oe 
shee os kp + ke’ Om 
(3-9) now yields 
cos 6, + cos 0, = 0, 
ee ee (3.28) 


cos 6, — cos 9, = 0, 
v=s,s+l,....,.st+tq-—l1. (3.29) 
Using (3-28) and (3.29), the last of the equations (3.25) becomes 
k, (cos? 6, + cos? 6, +....+ cos? 0) 
— ks (cos? 0, +....+ cos? 95,94) = 0, 
which together with the first of the equations (3.25) gives 
cos? 6, + cos? 8, +....+cos? 4, 
= cos? 8, + cos? 6, +....+cos?0, = ks/kp + ks, (3-30) 
cos? 8, + cos? 54, +....+ COS? Os49-4 
= cos? 8, + cos? O54, +....+ COS® Og4q-1 = Krlkp + ke. 
(3.24) and (3.30) show that the directions t and t are linear combinations of 


ee x €y, COS Oy, and A "S eo, cos 8,, kr - &s 
v=s $s 





the directions 
$s w=1 T 


r - ky a ks stq—1 ~ P ° . ° 

D eu, cos 6, and — x €y, COS Fy respectively, which being linear 
uel T o=s 
combinations of the principal directions corresponding to equal values of 
curvature, are themselves principal directions. 


Hence we find that the characteristic directions are linear combinations 
of the principal directions corresponding to any two distinct values of the 
principal curvatures. 


The angle between the characteristic directions t and t is given by 


ng — Am hoe (3.31) 


by virtue of (3.27). 
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The normal curvatures in the characteristic directions t and t are 
k = ky (cos? 6, +....+ cos? 6,) 
+ ks (cos? 0, +....-+- Cos* 9549-4) 
= 2krks/ky + ks. (3 -32) 
and 
k = ky (cos? 6, +....+ cos? 6,) 


+ ks (cos? 6, +....+ cos? Os1q-1) 
= Ukpkslkee + ks, (3-33) 
by virtue of (3.30) 
Relations (3 .31)-(3 .33) are the same as (3.20)-(3.22). 


4. DIRECTIONS FOR WHICH THE RATIO OF THE GEODESIC TORSION 
(oR SECOND CURVATURE) AND THE NORMAL CURVATURE IS EXTREMUM 


The geodesic torsion in the direction t= Y ep, cos 4p is given by 
(Manikarnikamma, 1952). z 


Tg" a 3 z z cos* Op cos? 9q (kp (% Kq)? (4-1) 
The normal curvature in the direction of t is 
kn = & kp cos? Op (4.2) 


The directions for which tg/ky is extremum are the same as those for tg?/k,y? 
is extremum, (for tg is zero only in a principal direction). 
Now, 

Tg*/ky? aaa (tg? + Ky?) ky? — 1 


= J k,? cos? 6, ei Kp cos? 9p) “2-1. (4.3) 
. Pp 
The directions for which tg/kp is extremum are therefore given by 


wer z k,? cos? 6, 


3 2p (5 ky cos? ®,), di 





1] = 0, 


v 


a 4.4 
2» cos? 6, = 1, os 


p=1 


and 


3 cos Op sin bp dy = 0. 


9=1 
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‘The first of the equations (4.4) may be written as 
= [kp? cos Op sin 0, (Lk; cos? 4,)? 
p 
—2 (Xk; cos* 6;) kp cos Op sin 0, 2'k,? cos? 9,] d8,=0 (4.5) 
Applying the Lagrange’s method of undetermined multipliers we have, 
( z k, cos? 6,) Ky* cos Op sin 8p — 2kp cos Oy sin 9 Sky? cos?6, 


+ Acos 9p sin 6, = 0, (4.6) 
(p = 1, 2,....n) 
where A is the Lagrange’s multiplier. 


Now sin 9, 40 for any p, for otherwise the direction t must be a principal 
direction. (4.6) therefore gives 


£08 Op | kp? 3 ky cos* 0 — kp 3 kr? cos* 8, + 2] = 0 (4.7) 


(p = 1, 2,....n) 
(4-7) gives either (i) Cos 6, = 0, 


or (ii) kp? 3 ky cos® 0, — 2k E ky? cos? 6, + A= 0. 


As in §3 two different cases arise: 
Case (i). When all the principal curvatures are unequal. 


Since t is not a principal direction, there must exist at least two values 
of p, say r and s, such that cos 6) = 0, i.e., for at least two values of p, 


k,? z Kp cos? 0, — 2kp B kp? cos? 6, ++A=0 (4.8) 
Putting p = r, s in (4.8) and subtracting we have 

(kp + ks) Eke cos? Op — 23 kp? cos? Oy = (4.9) 
since kr # ks. 


Multiplying (4.8) throughout by cos* 4, and summing up on p we have 
d= J kp? Cos? Oy 3 kp cos? Op (4.10) 
From (4.9) it follows that there cannot exist more than two values of p such 


that cos 8, may be different from zero; for, if r,s, w be three such values, 
we have as in (4.9) 
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(kp + kw) E kp COS? 0 — 2 3 kp? cos? Op = 0 
: 4.1 
(ks + Kw) 3 kp cos? 0, — 2 Y kp? cos* 0, = 0 
By subtracting the equations in (4.11) we find that either k, = ks or else 
z Ky cos* 0, = 0, 
ie., direction of t must be asymptotic. Since none of these alternatives is 


possible, we find that 
cos 6, = 0, p= 1,2,....4,p #7, 58, (4.12) 


ky? 3 kp COs? Op — 2kp 3 kp? cos® Op + A=0 
a A (4.13) 
ks? 3 Kp cos? Op — 2ks Y kp? cos? 6, + A=0 


The second of the equations (4.4) now gives 
cos? 6, + cos? 6, = 1 (4.14) 
Using (4.12), (4-9) gives 
tan? 0, = k;/ks, 
tan? 0, = ks/ky. j 


(4.12), (4.15) and (4.16) define a pair of real directions corresponding to 
every pair of distinct values of principal curvatures of the same sign. 


(4.15) 


CAsE (ii). When there are sets of equal values of principal curvatures at 
any point. 

Let ky = ky =....= kp Fikes = Kegyy =... . = Kgaq, ete. 
Proceeding as in § 3 above, we can show that the directions for which the 
ratio of geodesic torsion to normal curvature is extremum are given by 

cos 6) = 0, p #1, 2,....7,5,5+1,....8s+q—1, (4.16) 
k. 
2 ee... see 
cos? 0, + cos? 4, +.... + cos? 4, E+k 


k 
cos? O5+c0S? Os41+.... C08" Osiqa= F + Ke 


(4.17) 


Comparing (4.12), (4-15), (4.16) and (4.17) with (3.15), (3.19), (3-24) and 
(3.30) respectively we find that the directions for which the ratio of the geodesic 
torsion and normal curvature is an extremum are characteristic directions, 
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and are, as such, linear combinations of the principal directions corresponding 
to two distinct values of the principal curvatures. 


5. PARTICULAR CASE 


If the enveloping space be an Euclidean 3-space, we find from §4 above 
the known result (Mishra, 1947), that the characteristic lines on a surface 


are the directions for which the ratio of geodesic torsion and the normal curva- 
ture is an extremum. 


6. SUMMARY 


Upon a surface of positive Gaussian curvature there exists a unique 
conjugate system for which the angle between the directions at any point 
is the minimum angle between the conjugate directions at that point. This 
system of lines is called characteristic lines. In the present paper charac- 
teristic lines of a hypersurface V, imbedded in a Riemannian V»,, have been 
studied. It has been proved that characteristic directions are linear combi- 
nations of the principal directions corresponding to any two distinct values 
of the principal curvatures. It has also been proved that the normal curva- 
tures in the two characteristic directions lying in the pencil determined by 
the principal directions corresponding to two distinct values of the principal 
curvatures are equal, each being equal to the harmonic mean between the 
principal curvatures. The directions for which the ratio of the geodesic 
torsion and normal curvature is an extremum have also been studied and it 
has been shown that the directions for which the ratio of the geodesic torsion 
and the normal curvatures is an extremum are characteristic directions. 
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SPECTRUM OF AXI-SYMMETRIC TURBULENCE 
IN A CONTRACTING STREAM 


By Dr. Y. V. G. ACHARYA* 


Received April 16, 1956 


(Communicated by Dr. P. Nilakantan, F.A.sc.) 


1. INTRODUCTION 


BATCHELOR! considered the passage of isotropic turbulence in a ‘ sudden’ 
contraction. The condition ‘sudden’ ensures that the problem is a linear 
one and tractable. 


Further extensions were made in solving this problem by Ribner and 
Tucker.2 They considered the spectra and the turbulence levels at points 
upstream and downstream of a contraction and derived the relations between 
the various quantities. Thus, if the initial spectra and turbulence levels in a 
sudden contraction are known and if the conditions are isotropic, the final 
spectral forms as well as the turbulence levels can be determined. If it is 
assumed that the contraction is axi-symmetric and the mean velocity ratio 
is 1,, we have from considerations of continuity 


ol, lJ, = 1, (1.1) 
where o = 1, for incompressible flows and /, = /, for axi-symmetry. Ribner 
and Tucker have determined the form of the spectral tensor at a point B 
downstream of the contraction in terms of the spectral tensor at point A 
upstream, assuming that the decay is negligible. Thus, if I°,,*(«) is the 
spectrum tensor at the point A, where « is the wave number vector, 





Pag B(x) = a » [Fast ( x) + (ake ae Tyg) Ky sd see: Jy 


€k,? + Ko? + Ks” 
Tyy4xy?kakg (1 — €)? | 
(ex? + Kg* + x3")? 
Here, « = /,?//,2 and « <1 for large /, and 
Ky, Kay Kg = 1,Ky, 1,K3, 1,Ks. (1.3) 
(1.2) is the basic equation for us. This equation has been used by Ribner and 


Tucker, to obtain many results, assuming initial isotropy. In the following 
lines, we shall use the same equation to solve the more difficult problem of 





(1.2) 





* Senior Scientific Officer, Directorate of Technical Development and Production (Air), Minis- 
try of Defence. 
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initial axi-symmetry. In a wind tunnel, or other apparatus, where most of the 
aeronautical experiments are done, the geometry would suggest that, by and 


large, one encounters axi-symmetric turbulence rather than isotropic 
turbulence. 


2. THE INITIAL FORM OF THE AXI-SYMMETRIC TENSOR 
The general form of axi-symmetric tensor has been given by Batchelor, 
as 
$ij (Kk, A) = Axixs + BAA; + Cdij3 + Dejdj + Exjaj, (2.1) 
where A, B, C, D and E are the defining scalars and functions of «? and KA. 
For the spectral tensor I'yg («, A): 


Katagp—K* Tug =0 
B = Kp! ap } (2.2) 
Dag (x, A) = Iga (x; A) 
Following Proudman and Reid,*® we define the operators 
Aoi = 84: — “83's Agi = 89; — “BQ). (2.3) 


Then the general axi-symmetric tensor which is solenoidal in both indices 
a and 8 is given by 


Pag (x, A) = Aq! (x) A,’ () Piz («, A) 
=B (A en —) ( hy = “epi 


Kqkj Kgkj 
+ C (845 —“s5!) (85 -“8), 2-4 
where B and C are the defining scalars, which are functions of x? and «;Aj. 
Instead of B and C we write I’, and I’, for the defining scalars and rewrite 
the general spectrum tensor in the form: 
Drop (0%) = Ts (845 — “2e?) (895 — “B) + 


K2 
+1, (a — “*#"*) (a, — 299). 2.49 


Equation (2.4 a) has to be used in conjunction with equation (1.2) to obtain 
the form of the spectral tensor downstream of the contraction in our case. 


The evaluation of the integrals becomes very difficult as the defining 
scalars are functions of x? and «;A; in this case. Whereas in the case of 
isotropic turbulence, which was dealt by the preceding workers, the defining 
scalar is simple, as there is no preferred direction. 
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The diagonal terms of the tensor, given by (2-4 a), at the initial section 
A are: 


Ty (x, 4) = r(i- )+r(1-% ' 





Pad, )=T,(1-F)+ 7% oy (2.5) 
Pat (,) =0,(1-8 **) > EE, 


The other three terms of the tensor can also be evaluated and will be men- 
tioned later. 


3. SPECTRAL FORMS DOWNSTREAM OF THE CONTRACTION 


The form of the spectral tensor at the section B could be obtained by 
introducing the corresponding equations in (1.2), for the terms at the sec- 
tion A. Thus, 

T° (K, ) = 


Kt 


L 
4 (x) (ex,? + x? + Kg?) 


a Ty" 





(3.1 a) 


where I’,,“ [we will hereafter write I",, instead of I ap (x, A) as there will 
be no confusion] is given by the first line of equation (2.5). Hence, we obtain 


le? [r.(1 —S)+0, (i-* )'] a 





ry,® (K, A) = L (€x,? + Ke" + wr . (3 -l) 
Similarly, we get , 
Pig se 
Pa? (K,) = h [Pat + Dera — 9 4 


DPyyA«y?x? (1 — €)? 
(ex? + at eb 
where I,,“ and Ij.“ are given in equation (2.5) and 
Pag! (x, 9) = — “2 [7 +r, (i-< om :)}. (3.3) 
Hence equation (3.2), similar to (3.1), can be obtained by substituting the 
values of I°;*, I'y,*, I"sg* from equations (2.5) and (3.3) in (3.2 a). 


We can, in a similar fashion, obtain the value of I’;38 (K, A) from (3.2 a) 
by substituting Iy*, I';3“ and «gs instead of I’,,4, I',s4 and «.; also it can be 
derived directly from (1.2) and (2.4). 


(3.2 a) 
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4. DETERMINATION OF THE CORRELATION TENSOR 


The relation, between the correlation and the spectrum tensors, can be 
expressed by the following equations: 


Rast (7) = STS Top’ () ete? de (1) 


‘ei (4.1) 
Rag? (”) = JJ J ag? (K) e-* dr (K). 
The auto-correlation at the point is further given by 
+00 
Rea® 0) = SSS Pua (x) dr (1) 
(4.2) 


Rac? (0) = J Ff Poa? (K) dr (K). 


So, we have got a method of determining the auto-correlations of the 
turbulent fluctuations at the point B downstream of the contraction, if we 
know the spectrum at a point A upstream of the contraction. 


5. EVALUATION of u,? 


The turbulent fluctuations at a point are designated u, v and w as usual, 
then we have the following relationships: 


uy? = Ry (0); v4? = Reo! (0); wa? = Rgs* (0) (5-1) 
Uy? = Ry? (0); vg? = Ro»* (0); Wa? = R,j,° (0). (5.2) 
For convenience, in integration, we use spherical polar co-ordinates: 
kK, =«cos@ 
Ke = «sin 6 cos ¢ 
Ks = «sin 6 sin ¢ (5.3) 
dr(x) = x* sin 6 dé dé dk 
2 
dr (K) = 7 sin 0 d0 df de. 
hhls 


We shall now try to evaluate the individual mean squares of the turbulent 
fluctuations. 


(a) Variation of u® 
We have from (5.1) and (4.2), for the section A: 


u,? = Sis T'y;4 (x) dr (x), 
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where 
T',y4 («, 4) = I’, sin? 6 + I, sin @ (5.44) 
from the first line of (2.5), by introducing the values given by (5-3). 


Similarly, we obtain from (5.2) and (4.2), the value of uy? as 


— +co 
up = J J J Pu® (K) dr (K), (5.5) 
where 
__ 1,? (1, sin? 6 + I, sin* 4) 
Tu® (K, %) = Ta “(ecost 8+ sin® 6) * wii 
Hence, we have 
+00 

=a SJ J Pu® (K) dr (K) 
Sw a (5.6) 





ust FFP Tuk (e) de (x) 
(b) Variation of v® 

From the second equation (5.1) and the first line of (4.2), we obtain the 
value of v,?, as: 
mel +co 
DE = STI Tae! («) de (1), (5-7) 


where 
Tyo4 (x, A=", (1—sin? 6 cos? ¢)+T, sin? @ cos? @cos? ¢. (5.7 a) 


For the section B, vz” is obtained by using equations (5.2) and (4.2), as 











sy 
vp? = sJJ T'z2° (K) dz (K), (5.8) 
where 
Tn® (K, d) = |, [7 (1—sin? @ cos? ¢)-++I’, sin? 0 cos? 8 cos? ¢ — 
__ 2 sin? 6 cos? 8 cos? ¢ (lI—«) (71+ I; sin? 6) 
e cos? 6 + sin? 6 + 
4 sin? 8 cos? 6 cos? ¢ (1— «)? (I; sin? 6+, sin* 4) | 
(€ cos? 6 +- sin? 6)? 
5.84 
Hence, ¢ ‘ 
+00 
ae | § SJ Txg® (K) dr (K) 
== = (5.9) 


© 


2” iiroere 











68 Y. V. G. ACHARYA 





(c) Variation of w? 
The value of w? is obtained at sections A and B in a manner similar to 
ve 


the calculation of v*, for this purpose we have to use I; in place of I... 


Hence 


= Sf Ta®(K)dr(K) 
OSES Tas () dr (0) 


= 





| 
| 


(5-10) 


= 


6. EVALUATION OF THE INTEGRALS IN EXPRESSIONS FOR 4,2 


The general problem of the evaluation of the integrals (5.4) — (5.10), 
is very difficult as the defining scalars I’, and I, are functions of «? and «Aj, 
Also the forms of these functions are not known. But we can evaluate those 
integrals under the initial assumption that I, and I, are only functions of «, 
and obtain a specialised picture of the flow. In the following lines, we shall 
discuss the results under this assumption. 


(a) Variation of u* 


The integral expression (5.4), can be solved with the aid of (5.44), 
resulting in 


a = (®, + : ®,), (6.1) 
where, we have put 
©, =f Ty «de; 0, =f Ty «*de. (6.1 a) 
When, (5-5 a) is —_ in (5-5) the integral becomes 
~ ia [O,1, + ,1,), (6.2) 


where, 





l—e 

2+e«  €(4—86) 
aq—" a=" 
Hence, we get for (5-6): 


uz? _—sid33s FL, + OI (6.3) 
A 1+ se 


I,” =| =. 4 (=o arc tanh vi=*]) (6.2) 


an 
I, = 
























(0 


fr 
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(b) Variation of v 
Expression (5.7) can be evaluated by introducing in it the value of I’..“ 
from (5-7 a), leading to 


8 1 
v2 = 7 (% + 10 ,), (6.4) 


Similarly, equation (5.5) can be solved and simplified to: 





vp? = rs (PL, + %,1,], (6-5) 
2 
where 
ae 2 
1?) = : — = —a _ aye are tanh 1 —« 


— 2é e?(2 + ©) 
=" (= 9 


Hence, equation (5.9) becomes: 


=_. 3 be + oh 
de 242 | (6.6) 
®,+ 19% 


L{* = 


arctanh YW1l—e« (8-5a) 





(c) Variation of w? 


w? varies in the same way as v v® because of axi- symmetry and we can write 
the expression (5.10) as: 


wa 3 Oh” + OL 
= =— a7 13 242 (6.7) 
Wa ®, + i ®, 





where the quantities I,‘*) and I,*) are given in the same way as in (6.5 a). 


We would like to emphasize the fact that expressions (6.3), (6.6) and 
(6.7) are true only if I, and I, are only functions of «, but they would not 
hold, if they are functions of «? and «;A; in a general way. There is no way 
of predicting this, unless experiments are made and the spectra evaluated 
in the case of axi-symmetric turbulence. 


7. ONE-DIMENSIONAL LONGITUDINAL SPECTRUM 


If we make a Fourier analysis of the variation of velocity along a line, 
the spectrum functions can be obtained by integrating over all « for the 


lateral components. Thus, the contribution to the mean square velocity 
A2 
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component u,? [which is equal to R,,'], in the x,-direction, from all the 
wave numbers between | «, | and | x, | + | dk,| is 


@, (K,) dk, = [255 Pag dks dr | dk, 


or 
@, (1) = 25 f Pag dg dg (7.1) 


which yields 9,4 («,) or @,8(«,) according as we use I’,,“ or I’,,°, in place 
of I,. in the integration of the R.H.S. 


@, is called the one-dimensional longitudinal spectrum of the flow 
turbulence. Similarly, we can obtain the lateral spectra when a = 2, 3 in 


terms of the scalar functions I, and I, which define the axi-symmetric 
turbulence tensor. 


8. CONCLUDING REMARKS 


Just as we carried out the calculations in paragraph 4, we can also deter- 
mine the one-dimensional spectral functions, under the assumption that I’, and 
I, are functions of only «. But, since no purpose is served unless we are 
sure of the functional dependence of I’, and I’, on «? and «;A;, we do not want 
to take the trouble to solve (7.1). But the problem has been stated and 
once the forms of the defining scalars are determined, by experiment* or 
otherwise, the calculation becomes straightforward, though laborious. 


It is intended to tackle the non-linear problem of a contraction, which 
is not sudden, in a subsequent paper. 


9. SUMMARY 


The passage of axi-symmetric turbulence in a suddenly contracting stream 
has been investigated. A form of axi-symmetric tensor suitable for this 
study has been derived. Assuming that the form of turbulence is initially 
axi-symmetric, upstream of the contraction, the spectral forms downstream 
have been given. The ratio of the turbulence levels, on either side of the 
contraction, has been evaluated in integral form. Further, the integrals 
have been solved, assuming that the defining scalars are functions only of 
the wave number x. The method of evaluating one-dimensional spectra 
has been indicated. The whole investigation is based on the linearised 
approach. 


*Note: Dr. Mahinder S.Uberoi of the University of Michigan, U.S.A., wrote to the author, 
while the latter was at the Johns Hopkins University in July 1955, that he was carrying on 
some experiments relative to the problem treated in this paper. 




































Spectrum of Axi-Symmetric Turbulence in a Contracting Stream 71 


e 10. NOMENCLATURE 

u turbulent fluctuation. 

co-ordinates. 

h mean velocity ratio. 

kh vertical and lateral contraction ratios. 
) K wave number vector. 

| 0 compressibility factor (o = 1, for incompressible flows). 

, : r spectral tensor. 

R correlation tensor. 
; A,B sections upstream and downstream, respectively, of the 
" contraction. 

1,2,3; a, 8, y subscripts for denoting directions. 

r,, I; defining scalars. 
| e=1,%/l,?. 
; Other symbols are defined in the text, as and when they appear for the 
. first time. 
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ABSTRACT 


The overlap factors for bond electrons in various diatomic halides are 
calculated making use of a simple relation where the overlap factor is an 
inverse exponential function of the difference between the two atomic radii 
of the molecule concerned. These overlap factors together with the 
hybridization values of 15 percent for Cl, and 10 percent for Br and I 
are used to calculate the ionic characters for various diatomic molecules 
from the values of the nuclear quadrupole coupling constants obtained 
by earlier workers. These ionic characters are then plotted against 
electronegativity differences and the points are found to be fairly well 
represented by a smooth curve. Jt is shown that this curve can in general 
be used for the determination of the ionic characters in the halogen con- 
taining single bonds of certain polyatomic molecules. The quadrupole 
coupling constants are calculated from the ionic characters obtained from 
the curve and are compared with those observed. The effects of the dis- 
tortion of the closed shells and of the neighbouring ions on eQgq are 
discussed. 

INTRODUCTION 


THE purpose of this paper is to discuss the method of obtaining the ionic 
character in halide molecules in gas phase from the nuclear quadrupole 
coupling constants.!>? These coupling constants which are obtained from 
the quadrupole hyperfine structure are denoted by eQg, where e represents 
2 
the electronic charge, Q is the nuclear quadrupole moment and q = ae 
is the gradient of the electric field at the nucleus due to extranuclear charges 
and depends primarly on the electrons in the p orbit of the valence shell; 


z is the co-ordinate along the symmetry or bond axis. 

The quadrupole coupling constant for a halogen atom forming a hybrid 
bond with a negative ionic character of amount 8, is given by® 

CQFmoieeure = (— 1 + 5 — d — S*) (1 — B) CQG prom (1) 

where s and d represent the respective contributions of the corresponding 
electrons in the bonding orbital of the halogen atom which is given by 

* Aligarh Muslim University Research Scholar. 
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Prtstocen = (1 — 5s — dt Pp + st Ps + dt Pq 


§ represents the overlap factor |¥,'%,"d7 for two atoms A and B of the di- 
atomic molecule, and eQq,,,,, is obtainable from atomic data.’ 


Since it is very hard to distinguish quantitatively the effects of s character 
from those due to d character, it is not necessary to retain separate symbols 
for them as in (1), and one may use the symbol a for the combination (s — d). 


If the bond is positively ionic that is when the halogen ion is positive, 
two electrons are missing from the valence shell of the halogen atom, and 
therefore the number of unbalanced p electrons will increase to 2, giving 
rise to a gradient® of the electric field equal to 2q+. Here g* represents the 
gradient of the electric field due to one valence p electron when the atom. is 
positively ionic, and is somewhat greater than the corresponding ‘ q’ when 
the atom is neutral. This is because the interaction between a valence p 
electron and the nucleus is somewhat larger in this case than in the case of 
the neutral atom. The ratio of q+ and q is taken as (1 + «), where « is nearly 
0-25 for most atoms.* Hence the eQq of a halide molecule for a positive 
ionic state will be given by — 2(1 + €) eQqa,,n, Where the hybridization has 
been neglected because of the increase of promotional energy required. 
Therefore, the quadrupole coupling constant for an atom with a positive 
ionic character of amount f is given by 


CQ4nor. = (—1 + a — S*) (1 — B) — 21 + 2) B) CQGatom (2) 
The exact value of 6 can be obtained from (1) and (2) if the values of a and 


S* are known. Gordy®? assumes that both a and S? are negligible, but the 
following difficulties arise with such an assumption :— 


(1) The ionic characters for FCl and FBr with such an assumption come 
out to be about 26 percent and 33 percent respectively* from equation (2), 
whereas from Gordy’s graph’ they should be about 50 percent and 60 per- 
cent respectively. Gordy says that this inconsistency could be avoided if 
one assumes hybridization in Cl and Br atoms of FCl and FBr in the positive 
ionic state. To explain this inconsistency one has to assume the hybridiza- 
tion for positive ionic importance to be 31 percent and 28 percent in FCI 
and FBr respectively,f as the hybridization has to be taken as negligible for 
covalent importance according to his views. This seems to be improbable. 

* From Gordy’s formula (21) (ref. 7) the values of ionic character come out to be 28 percent 
and 35 percent for FC] and FBr respectively. 

+ If Gordy’s assumption is correct, formula (2) has to be rewritten as 

€Q4 mor. = — [A — A) +20 + € — a) B] CQGg ton 


If g for FCl and FBr has to come out as 50 percent and 60 percent respectively from 
this formula, a has to be taken as 31 percent. and 28 percent respectively. 
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(2) The ionic character for BrCl molecule considering the coupling 
constant for Cl nucleus comes out to be 5-6 percent, whereas it comes out 
to be 12 percent while considering Br nucleus. 


(3) Gordy’ points out that in obtaining ionic characters for diatomic 
molecules from dipole moment data, one need not find it necessary to take 
hybridization into account, which he says is in support of his assumption that 
hybridization is negligible. But one should remember that the hybridiza. 
tion effect on nuclear quadrupole coupling constant is different from that on 
molecular dipole moment. Moreover, the effects of hybridization of the 


two atoms in a diatomic molecule act in opposite directions for dipole moment 
considerations. 


Dailey and Townes? took a into consideration and neglected S?, because 
no satisfactory calculations of the overlap effects seemed feasible. They 
followed a simple rule that the halogen bonds have 15 percent s character 
whenever the halogen is bonded to an atom which is more electropositive 
than the halogen by as much as 0-25 unit, and 0 percent otherwise. If one 
generalizes this rule and makes it applicable to Cl in BrCl where Br is more 
electropositive than Cl by 0-2 units, the value of 8 neglecting S* comes out 
to be negative which does not have a physical meaning. It is probably for 
this reason that they restricted 15 percent hybridization to cases where the 


halogen under consideration is more electronegative than the bonding atom 
by more than 0-25 unit. 


The hybridization in a positively ionic state is suppressed due to in- 
crease of promotional energy and therefore it may be assumed to be 0 per- 
cent for the positively ionic importance only. But Dailey and Townes* 
assumed 0 percent hybridization even for covalent importancet in their 
calculations on FCI, FBr, etc., all of which have only partial ionic character, 


In this present paper both a and S? will be taken into consideration. 
The values of a and S? will be based on certain general assumptions keeping 
the physical ideas in mind, as it would be difficult to obtain these values more 
precisely. In this paper the assumptions of Dailey and Townes® regarding 
hybridization will be modified, so that the difficulties pointed out above will 


not arise. Further the modified assumptions could also be extended to the 
case of polyatomic molecules. 


{ This they do giving the reason that if the resonance energy between the covalent and ionic 
states is to be maximum, it is necessary that both states possess the same amount of hybridization. 
If this reason is correct the hybridization has to be taken as something between 0 and 15 percent 
and certainly not 0 percent because in the molecule considered, the covalent character is more 
predominant than the positive ionic character. 
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ASSUMPTIONS REGARDING OVERLAPPING AND HYBRIDIZATION 


(a) Overlapping in diatomic molecules:—Overlapping is essentially a 
measure of the amount of interpenetration of the bond-electron distribu- 
tions of the two atoms,® and has a finite value whenever a bond exists between 
the two atoms. Thus one can start with the assumption that the overlap 
factor S is a decreasing function of either Rap or | ra — rp |, where Rap 
represents the distance between the two nuclei, and rq and rp represent the 
radii of the two atoms forming a diatomic molecule. Now using the atomic 
radii given by Slater?® and assuming Rap = (ra + rp), the overlap factors for 
various homopolar molecules are obtained from the tables given by 
Mulliken et al.,1° and are given in Table I. Here only those molecules in 
which the atoms have one electron in the outermost orbit of ground state 


TABLE I 


Overlap Factors in Homopolar Molecules 








Molecule Overlap factor R,, nA |ro—T,|inA 
i. & me 0-407 2:24 0 
x Sue 0-385 1-92 0 
3. Al, i 0-385 2:74 0 
4. Na, te 0-375 4-34 0 
> Ff, + 0-320 0-82 0 
. a 0-320 1-64 0 
% iy re 0-456 3-26 0 
8 Hs eh (0-586) 1-06 0 





were considered and it was assumed that the overlapping takes place only 
between these two electrons which are also assumed to form the covalent 
bond. 


This table (I) shows that the variations of S with Rgp are small, and the 
values in different homopolar molecules are fairly constant. Thus it appears 
that S can be more correctly represented by a function of | rg — rp | than 
by that of Rap. Therefore assuming that S becomes negligible as | rg — rp | 
becomes extremely large, a simple functional form of S? is 


3 Ae-2 (IT .-T,!) (3) 
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where A is a constant. The value of this parameter A is found to be 0-2] 
from the coupling data for BrCl, taking the ionic character f to be 6 per- 
cent," and assuming a for Cl in this molecule to be 15 percent.*» 12 This 
value of A introduced in (3) yields S = 0-45 for a homopolar molecule in 
fairly good agreement with the average value of S from Table I. 


(b) Hybridization —When the bond is purely covalent, a for the halo- 
gen atom in question arises due to the mixing of the valence p orbital with 
the s or d orbitals of the same atom. It will be fairly correct to assume that 
for the covalent importance the a for a particular halogen atom has the same 
value in various molecules of the same halogen, but may have different values 
for different halogens; for example, a for Br and I will probably be smaller 
than for Cl, because the d character in a for Br and I is greater than that 
for Cl. 


The promotional energy required for positively ionic state of the atom 
will be increased and hence the hybridization is reduced, so that it can be 
assumed that for this positively ionic importance the hybridization is sup- 
pressed to 0. 


Considering Cl in IC] and assuming 15 percent a for Cl, the ionic 
character in the molecule was found to be 21 percent. | Now considering 
the iodine atom in IC] and keeping § as 21 percent the value of a for I is 
found to be 9 percent. An approximate round figure of 10 percent will 
therefore be taken as the value of a for I in different molecules. Following 
a similar procedure the value of a for Br in different molecules is taken as 
10 percent. 


Incidentally it may be mentioned that the above assumptions would 
enable one to predict the eQq in the molecules I,, Br, and Cl, in gaseous 
phase to be — 2420 mc., 811 mc. and — 109-5 mc. respectively. * 


CALCULATIONS OF IONIC CHARACTERS IN DIATOMIC HALIDES 


The values of S? have been calculated for a number of diatomic molecules 
from (3) using the atomic radii given by Slater,® and are included in Table II. 
Using these values of S*, and assuming a for Cl to be 15 percent and that for 


* These values of eQq for these homopolar molecules are calculated from the equation 


Qua), = — I — a) (2 — 7g) Qtaren 


which, when S?is small, reduces to the equations (1) and (2), where g is to be taken as 0 for homo- 
polar molecules. 


The eQg in gas phase for these molecules are not observed so far but the corresponding 
observed values? in solid state are — 2153 mc., 765 mc. and — 108-95 me. for I,, Br, and Cl, 
respectively. 
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Br and I to be 10 percent, the ionic characters for the various molecules 
have been calculated from (1) and (2), and are also given in Table II. 
These values of the ionic character obtained are plotted against the electro- 
negativity differences (| xq — Xp |), and the shape of the resulting curve is 
shown in Fig. 1. The values of (| xa — Xp |) used for all except TICI are 
those given by Huggins.’* Since Huggins has not given the value of electro- 
negativity for Tl, the corresponding value calculated by Gordy® from a force 
constant relation has been used. 
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Tonic CHARACTER IN ALKALI HALIDES 


From equation (1) one expects that the magnitude of the quadrupole 
coupling constant for halides molecules should decrease asymptotically to- 
wards zero with the increase of ionic character which in turn is an increasing 
function of | xq — Xp |. On the other hand as the electronegativity difference 
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increases the quadrupole coupling constant in the chloride molecules, as 
shown in Table II, decreases in magnitude up to KCI where it changes its 
sign and increases in magnitude from KCl to CsCl. This change in alkali 
chlorides can be explained in the following manner :-— 


The contributions to q in molecules come from, 
(a) the valence electrons associated with the atom under consideration, 
(5) the neighbouring ions and atoms, and 


(c) the distortions of non-bonding closed shells of electrons around 
the nucleus in question, 


of which the first is the most predominant. If D represents the effect on 
quadrupole coupling constant due to distortions of the closed shells around 
the halogen nucleus, and i represents the effect due to presence of adjacent 
ions or atoms, the equations (1) and (2) for coupling constants can be 
more precisely rewritten as 


CQ4 nai. = [(— 1 + a — S*) (1 — B) + (D — i) B) CQGa tom 


(negatively ionic state) (4) 
CQ4mn = (— 1 + a — S*) (1 — B) —2(1 + &) B — (D — i) B] eQy,,,,, 
(positively ionic state) (5) 


From (4) it is clear that when 8 approaches 1, the term (D— i) 8 
€Q4 atom becomes appreciable and will be responsible for the change in sign 
of the quadrupole coupling constant if the distortion effects are greater than 
those of adjacent ions. This explains the change in sign of eQg for KCI 
and its increase in magnitude from then on up to CsCl. The best values of 
8 as obtained from the curve in Fig. 1 are 99-5 percent for KCl and RbCl 
and 100 percent for CsCl. With these values of 8 the factor (D — i) B 
CQAarom CaN be calculated to be equal to 0-44 mc. for KCl, 1-24me. for 
RbCl and 3-00 mc. for CsCl. As 8 ~1 here, this shows that (D — i) is 
increasing with the increase in interionic distance, which is quite under- 
standable as the term i decreases rapidly with the increase in interionic dis- 
tance. The presence of the factor (D — i) in (4) explains why quadrupole 
coupling constants as calculated from the simple formula (1) for alkali 
halides come out to be higher in general than those observed. 


It is interesting to see that, in the case of FBr where Br is partially posi- 
tively ionic and where the ionic character is about 40 percent the quadrupole 
coupling constant as calculated from (2) comes out to be also higher than 
that observed. This is understandable from (5) which is applicable to FBr 


if the factor (D — i) is taken as negative. This would mean that the effect 





TABLE II 
Ionic Characters and eQq for Diatomic Halides in Gaseous Phase 
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of the neighbouring ion in the ionic state of the molecule is larger than the 
distortion effect, whereas it is smaller in alkali halides, which is reasonable 
because the interionic distance in FBr is much smaller than that in any 
alkali halide. It is to be pointed out that the effect of the term (D — i) ig 
also noticeable in FCI in a similar way as in FBr but in a smaller degree as 
8 for the former is smaller than for the latter. From the deviations between 
the observed and calculated values of eQg, it can be estimated that the effect 
of the term (D — i) in the ionic states of FCI and FBr is to decrease the 
quadrupole coupling constant by about 8 percent and 12 percent respec- 
tively. Thus the combined effect of neighbouring ions and of the distortion 
of closed shells is in general to decrease the eQg in diatomic molecules, 


IONIC CHARACTER IN POLYATOMIC MOLECULES OF HALOGENS 


The ionic characters for the single bonds containing halogens in various 
polyatomic molecules have been calculated from (1) neglecting the overlap 
effect} and the corresponding values obtained are given in Table III. These 
values are plotted in Fig. 1 and are marked with cross (x). It is interesting 
to see that all the points in general fit in fairly well with the curve drawn for 
the diatomic halides. 


In the case of polyatomic molecules, the effect of the involved double 
bond character! 1° on eQq should be taken into consideration besides the 
other effects discussed for diatomic halides. In this paper we are, however, 
considering molecules involving a small percentage of double bond character 
only. The effect of the double bond character is to increase eQgq,'* but the 
combined effect of the distortion of closed shells and the neighbouring ions 
as mentioned earlier is to decrease eQqg. The fact that the ionic characters 
calculated from (1) fit in general with the curve in Fig. 1 shows that the 
effect of the double bond character is approximately compensated by that of 
the term (D — i). Therefore, it appears that the ionic character for a single 


+ Considering a molecule like H,C—Cl one can see that because the carbon atom for co- 
valent importance forms single bonds with each of the other four atoms of the molecule, the atomic 
orbital , of Cl will have certain possibilities of overlap with each of the bonding orbitals 


Pos Pos. and Pon, of carbon. These possible overlap factors are proportional to 
|.» Psy, dt, [Pep tn! "tg dz and JP iy P29, dt. 







But 

JP as Pyg dt 1/3 [Pang Pan, de® and [Pop 4g Payg dt 
can be shown to be zero. Thus the average value of the overlap integral for the C-Cl bond is 
small compared with f¥,,, ¥*;,, dz which usually represents the overlap factor for a diatomic 


molecule. S? may therefore be neglected for approximate purposes in such kind of polyatomic 
molecules as considered in this paper. 
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The Ionic Character of Singly Bonded Molecules 


TABLE III 
Ionic Characters and eQq for Polyatomic Molecules in Gaseous Phase 














Ionic Ionic 
B Molecule Hybrid- Character Character /x,—-x,/ Observed Calculated Ref. 
ization calculated in % (Huggins) eQq eQq 
in % from graph in me. in me. 
h35 a is me 7 —109-74 =e oe 
,c-Cl 15 19-5 21 0-55 —75:3 —73-7 a,b 
SC] .. 15 57 54-5 he25 —40-0 —42°4 d 
Ge-Cl .. 15 51 54-5 1-25 —46-0 —42°4 g 
i & ie a ees a 769-76 ah 
C-Br 10 16-5 13 0-35 577°3 602-7 a,c 
i-Br 10 51°5 44 1-05 336-0 388-0 e 
Ge-Br .. 10 45 44 1-05 380-0 388-0 h 
i _ <2 Sy tg eS —2292-84 ot 
C-] 10 6 1-5 0-05 —1934 —2032 a 
ie 10 40 30 0-75 —1240 —1448 £ 
a Gordy, Simmons and Smith, Phys. Rev., 1948, 74, 243. 
6 Karplus R. and Sharbaugh, Jbid., 1949, 75, 889, erratum 1449. 
c Sharbaugh, A. H. and Mattern, J., Jbid., 1949, 75, 1102. 
d Sharbaugh, A. H., Jbid., 1948, 74, 1870. 
e Sharbaugh, Bragg, Madison and Thomas, Jbid., 1949, 76, 1419. 
f Sharbaugh, Heath Thomas and Sheridan, Nature, 1953, 171, 87. 
g Dailey, Mays and Townes, Phys. Rev., 1949, 76, 136. 
h Sharbaugh, Prichard, Thomas, Mays and Dailey, Jbid., 1950, 79, 189. 





bond involving a halogen atom in a polyatomic molecule can be approxi- 
mately obtained from the curve drawn for diatomic halides. The ionic 
characters obtained by this procedure are given in Table IV for the single 


TABLE IV 


Calculated Ionic Characters and eQq for Certain Polyatomic Molecules, 
the eQq for which is not so far observed in the Gaseous Phase 








Ionic eQq eQq 
Molecule Character in in me. in me. Ref. 

: percentage calculated for observed in 

from graph gas phase solid state 
SnCl, rae 54:5 —42°-4 —48-2 a 
PCI, bei 41 —55-0 —52°-3 a 
AsCl, ei ds —52°4 —50°-5 a 
SbCl, a 46 —50°4 —40-2 a 
Snl, Sa 30 —1448 —1363 b 





a Livingston, R., See ref. 1, p. 366. 
b Dehmelt, H. G., Naturwiss., 1950, 37, 398. 
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bonds in different polyatomic molecules along with the corresponding calcu- 
lated eQg. As the coupling constants for these molecules in the gaseous 
phase are not yet known, those observed for the molecules in solid phase are 
included in the table for comparison. 























The authors wish to express their thanks to Dr. A. N. Mitra for helpful 
discussions and to Professor P. S. Gill for his kind interest in the work. 


REFERENCES 
1. Gordy, Smith and Microwave Spectroscopy, John Wiley & Sons, Inc., 1953. 
Trambarulo 
2. Townes, C. H. and Ibid., McGraw-Hill Co., Inc., 1955. 
Schawlow, A. 
3. Dailey, B. P. and J. Chem. Phys., 1955, 23, 118. 
Townes, C. H. 
4. Jaccarino, King and See ref. 3. 
Stroke 
5. Townes, C. H. and J. Chem. Phys., 1949, 17, 782. 
Dailey, B. P. 
Gordy, W. .. J. Chem. Phys., 1951, 19, 792. 
Tt  cactemnieeaiilnaomes .. Faraday Society Discussion (Microwave and Radio-frequency 


Spectroscopy), April 1955. 


(The authors wish to express their thanks to Professor Walter Gordy for sending a manuscript 
copy of his discussion.) 


8. Pauling, L. -. Nature of Chemical Bond, Cornell University Press, Ithaca, 
1948, p. 76. 
9. Slater, J. C. .. Phys. Rev. 1930, 36, 57. 
10. Mulliken, Reike, Orloff J. Chem. Phys., 1949, 17, 1248. 
and Orloff 


11. Honig, Mandel, Stitch Phys. Rev., 1954, 96, 629. 
and Townes 
12. Paul N. Schaltz .. J. Chem. Phys., 1954, 22, 755. 
13. Huggins, M. J. .. J. Am. Chem. Soc., 1953, 75, 4123. 
14. Townes, C. H. and J. Chem. Phys., 1952, 20, 35. 
Dailey, B. P. 
15. Sheridan, J. and Ibid., 1951, 19, 965. 
Gordy, W. 
16. Mays, J. M. and Ibid., 1952, 20, 1695. 
Dailey, B. P. 




















COMPARATIVE STUDIES OF ALTERNATING 
AND DIRECT CURRENT POLAROGRAPHY : 
EFFECT OF pH 
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(Department of Physical Chemistry, Indian Institute of Sugar Technology, Kanpur) 


Received May 31, 1956 


(Communicated by K. S. G. Doss, F.A.sc.) 


ANALYSIS of a mixture of two or more electro-oxidisable or reducible sub- 
stances with their characteristic half-wave potentials, in a given supporting 
electrolyte, close to each other, cannot be carried out on the conventional 
direct current polarography; thus a mixture of Cd** and In**+ which have 
their half-wave potentials as —0-610 and —0-565 volts vs. S.C.E. could not 
be analysed on dc. polarography. This appeared to be a limitation to the 
polarographic work using dc. potentials. This difficulty could, however, be 
overcome by using alternating current polarography in which a known ac. 
voltage was superimposed over the dc. potentials fed to the dropping mercury 
electrode (d.m.e.) and the alternating current passing through the system was 
measured as a function of the applied dc. potential; these studies led to the 
observation of well-defined peaks due to the reversible reduction or oxidation, 
at d.m.e., of the substances, with their summits occurring at the potentials 
corresponding to the characteristic half-wave potentials recorded on dc. polaro- 
graphy and with the peak or summit current being proportional to the con- 
centrations of the electro-active substances." * A significant feature of these 
studies was that the alternating polarograms were remarkably sharp in the 
sense that the potential width of the waves was small, at low superimposed 
ac. voltages, say 5-10 millivolts (500 cycles/sec.), which led to the possibility 
of the analysis of the mixture containing Cd** and In*+++.1 Further, com- 
parative studies of the potential-widths of the dc. and ac. polarograms ob- 
tained at different concentrations of the electro-active substance also revealed 
that the alternating current polarography had remarkable advantage over the 
dc. polarography to the effect that by using low concentrations the former 
can be suitably employed for analysis of mixtures of two components with 
their half-wave potentials close to each other. These studies have now been 
extended to the investigation of the effect of pH on the wave-widths of ac. 
and de. polarograms of cadmium, on which no data exist in the literature. 
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EXPERIMENTAL 


The apparatus and the electric circuit employed have been described in 
the previous communication.® 


The dropping mercury electrode had the following characteristics: at 
a pressure of 32 cm. of mercury, drop time (¢) was 1-6 secs. in 0-1 M KCl 
solution with open circuit; and the weight of mercury dropping per sec. 
(m) was 2-5 mg. 

In the alternating current polarography, 30 millivolts (50 cycles/sec.) 
was used, throughout the present series of experiments, as the superimposed 
ac. voltage. 


Cadmium sulphate used was Merck sample and was recrystallised several 
times in twice distilled water and later dried at 160° C. 


The pH of the solution of Cd** of a fixed concentration was altered by 
using citrate buffers. Potassium citrate used for this purpose was of B.D.H. 
quality; hydrochloric acid and sodium hydroxide employed were Merck 
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Fic. 1. D.C. polarograms of Cd++ at different pH values. 
Curves 1 to 8 refer to pH values 4 to 11 respectively. 
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analysed samples. Required quantities of citrate and hydrochloric acid or 
sodium hydroxide solutions were mixed in the manner described by Meits* 
and others, with cadmium solution to give a sample of desired pH; the mix- 
tures were so prepared that the final solution was of the same volume and con- 
tained fixed quantities of Cd*+* and citrate. The pH of the solution was 
measured by a Beckmann pH meter. 


RESULTS AND DISCUSSION 


Figure 1 gives a typical series of dc. polarograms of Cd** at different 
pH values; these data refer to a concentration of 2 mM of Cd++ and of 0:5 M 
of Cit--- The corresponding ac. polarograms are returned in Figs. 2 (a) and 
2(b). It was interesting to note that at all pH values below 10, Cd** gave 
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VOLTS vs. S.C.E. 
Fic. 2(a). A.C. polarograms of Cd*+ at different pH values. 
Curves 1 to 4 refer to pH values 4 to 7 respectively. 
thermodynamically reversible waves on dc. polarography (Fig. 1) and well- 
defined peaks on ac. polarography (Fig. 2). At pH values greater than 10, 


however, no wave could be observed (Fig. 1) indicating that no reversible 
A3 
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VOLTS vs.S-C-E. 
Fic. 2 (6). A.C. polarograms of Cd*+ at different pH values. 
Curves 5 to 8 refer to pH values 8 to 11. 
reduction of Cd*+* occurred at the dropping mercury electrode‘; similar 
results were recorded on the alternating current polarography. This obser- 
vation suggested a convenient method for eliminating the interference of Cd** 


in alternating and direct current polarographic determination of other sub- 
stances.* 


Further, at pH values below 10, the diffusion current ig of the dc. polaro- 
grams, as well as the summit current i, of the ac. polarograms depended 
markedly upon the pH of the solution. Fig. 3 gives the variation of ig and is 
with pH. It was interesting to note that ig and i, decreased initially with 
increase of pH from 1-5. Thus ig was 26, 20 and 14 amps. at pH values 
1, 3 and 5 respectively (Curve 1, Fig. 3). At the pH values between 5-7, 
iq and is remained constant with pH. With further increase of pH above 8, 
both ig and i, decreased sharply (see Fig. 3). These features indicated that 
at varied ranges of pH, different characteristic equilibria leading to the exist- 
ence of thermodynamically reversible complexes of cadmium and citrate ions, 
were in operation. This was substantiated by the variation of half-wave 
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Fic. 3. Variation of diffusion and the summit current with pH of the solution. 

potential (E,) of the de. polarograms and the summit potential (E;) of the ac, 
polarograms, with the pH of the solution; these data are shown graphically 
inFig.4. It was interesting to note that at all pH values investigated, E, and 
E, were roughly the same as evidenced by the observation that the points 
representing E, and E, at varied pH fell roughly on the same curve (Fig. 4). 
Further, E, and E, increased with pH. In the region of pH values between 
5-7 the rate of increase of E, and Es was less than that recorded at pH values 
between 1-5; again at pH > 8, the rate of variation of E, and Es; with pH 
was considerably enhanced. The following processes appear to be occur- 
rent in the corresponding pH ranges: 


(1) Cd + Cit- - -=Cd Cit- + 2e (pH < 5) 
(2) Cd+* + H Cit.- - = (Complex) + Ht (pH ~ 5—7) 
(3) Cd Cit.- + OH- = Cd (OH) Cit. - (pH > 7) 


Enough evidence in support of these has been adduced by Meits,’ Lingane® 
and others*-*; the values of the corresponding equilibrium constants‘ have 
been evaluated ; and these are: K, = 6x 10°, K, = 5x 10+ and K, = 9x10-*. 
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Fic. 4. Effect of pH on the half-wave and summit potentials of Cdt+t. 


The data in Fig. 5 show the variation of the potential-width of the de. 
and ac. polarograms returned in Figs. 1 and 2. As indicated in the previous 
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communication® the wave-spread of the ac. peaks was taken as the potential- 
width of the polarograms corresponding to the current values equal to half 
the summit current. It was interesting to note that in the entire range of the 
pH values studied, the wave-spread of the ac. polarograms (Wc) was less 
than that corresponding to the dc. polarograms (W4q-). Thus, e.g., at pH = 2, 
the dc. polargorams spread over a potential region of 140 millivolts while the 
ac. waves occurred over 98 millivolts; similarly at pH = 10, Wg. and Wage 
were 85 and 60 millivolts respectively. The difference between the wave- 
widths of ac. and dc. polarograms was maximum at low pH values. These 
data indicate that analysis of a mixture containing Cd*++ as one of the com- 
ponents can be carried satisfactorily on ac. polarography in the pH range 5-6 
where the current values are appreciable and the wave-width is considerably 
smaller than in de. polarography. 


SUMMARY 


Alternating and direct current polarograms of cadmium were recorded 
at different pH values in the range 1-12. At pH < 11, cadmium was reduced 
reversibly at the dropping mercury electrode, giving well-defined waves on 
both ac. and de. polarography. The wave-width of ac. polarograms was less 
than that of dc. polarograms in the entire range of pH studied. 
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4-(p-NITROPHENYL-AZO)-8-HY DROXY-QUINOLINE* 
By N. A. RAMAIAH AND V. SUNDARA RAMAN 
(Department of Physical Chemistry, Indian Institute of Sugar Technology, Kanpur) 
Received June 6, 1956 
(Communicated by K. S. G. Doss, F.A.sc.) 


Azo dyes have been employed extensively as analytical reagents for metal 
ions. Gutzeit and Monnier’.* studied a large number of azo dyes containing 
the oxine nucleus and concluded that they gave specific tests for Hgt+, Pd+, 
Ni*, etc. Bogel Desering and Shreve* estimated the relative values of the 
azo dyes as analytical reagents for the varied metal ions. Of these azo dyes, 
4-(p-nitrophenyl-azo)-8-hydroxy-quinoline in alcoholic medium appeared to 
be important as an analytical reagent expecially for Ca*+.4> During an 
examination of its applicability to the estimation, with ease, of calcium content 
in sugarcane juices it appeared desirable to determine the dissociation con- 
stant(s) of the dye, the knowledge of which is of marked importance for clear 
understanding the use of the dye in analytical work. No data on this aspect 
exist in the literature. The present communication reports the determina- 
tion of the dissociation constant(s) by employing the characteristic absorption 
spectra of the substance in the visible region. 


EXPERIMENTAL 


The starting material employed by the authors for the preparation of the 
azo dye were acetanilide and 8-hydroxy-quinoline; and these were of the 
following quality: 8-hydroxy-quinoline used was of B.D.H. quality. Acetani- 
lide was a Bush product; and was converted into p-nitroaniline in the usual 
manner*; this last was purified by recrystallisation in rectified spirit till it 
gave a constant melting point 147-5° (148°) C. p-nitroaniline (1 mol.) was 
dissolved in N-HCI acid (24 mol.). The solution was cooled to 5° C. and then 
treated with 10% solution of NaNO, (1 mol.). The diazotised solution was 
added dropwise to alkaline solution of 8-hydroxy-quinoline (24 mol. in 
N-NaOH) maintained at 5° C.; the reaction mixture was stirred for one hour 
to ensure complete coupling.® The dye was then precipitated by acidifying 
it with 2 N acetic acid; filtered and washed several times with water. It 
was crystallised 4-5 times in rectified spirit; a product with a melting point of 
258° C. was obtained. 


*Apreliminary note was published in Naturwiss., 1956, 43, 179, 
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The absorption spectra of the alcoholic solution of the azo dye at different 
pH values was studied with a Unicam S.P. 350 D.G. Spectrophotometer ; 
this had one cm. effective light path. The buffers used were due to Clark 
and Lubs and to Sérensen’; the actual pH of the solution (1 vol. of solution 
+1 vol. of buffer) was checked by a Beckmann pH meter. 


RESULTS AND DISCUSSIONS 


Figure 1 gives a typical series of results on the absorption spectra of 4-(p- 
nitrophenyl-azo)-8-hydroxy-quinoline in alcoholic solution at different pH in 
the range 4-12. These data refer to a concentration of 0-002 mM of the azo 
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Fic. 1. Absorption spectra of 4-(p-nitrophenyl-azo)-8-hydroxy-quinoline. 
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dye. It was interesting to note that the dye exhibited absorption maximum at 
A = 600 mp in alkaline solution and another maximum at A = 480 my in the 
neutral and acidic solutions. It is to be remarked that freshly prepared soly- 
tion of the dye exhibited an absorption maximum at 560 my in alkaline solu- 
tions, which on long standing, however, shifted permanently to 600 mu. The 
experiments were conducted only when the system exhibited stable and per- 
fectly reproducible characteristics. 


It was interesting to observe that, while the absorption at A = 480 mp 
was not affected appreciably, the same at A = 600 mp was, however, variant 
markedly with the change of the pH of the solution especially in the range 5-12 
pH. The system exhibited an isobestic point, though not well defined in 
character, in the neighbourhood of A = 510 mp (Fig. 1). The data in Fig. 2 
show the variation of the optical density of the system at A = 600 my with pH 
of the solution. It was interesting to note that at A = 600 mp the optical 
density was almost unaltered with pH from 4-6-5; and a further change in 
pH from 6-5-7, increases the optical density markedly. Thus, e.g., optical 
density was 0-05, 0-08 and 0-2 at pH = 6-0, 6-5 and 7-5 respectively. Later 
in the pH range 7-5 and 8, the optical density of the system did not vary appre- 
ciably; in fact, it tended to saturate in the sense that the curve representing 
the variation of optical density with pH appeared to run parallel to the pH 
axis (see Fig. 2). Further increase of pH again increased optical density, 
which was unaffected by pH > 9-0 (see Fig. 2). Thus there appeared to 
exist two inflexions in the curve representing the variation of optical density 
with pH of the solution, indicating that 4-(p-nitrophenyl)-8-hydroxy-quinoline 
possess two dissociation constants. The data given in Fig. 2 are therefore not 
amenable for the familiar calculations, from spectrophotometric data, of 
the dissociation constant of the substance, which take into consideration that 
the optical density at different hydrogen concentrations is dependent on the 
proportion of the substance in the acid or base form (molecule or ion). The 
present circumstance is apparently complicated as in the case of the calcula- 
tion of the dissociation constant of a dibasic acid by potentiometric titration 
method.? The existence of two dissociation processes in the system under 


investigation is presumably due to the presence of the -OH and > N in the 
molecule of the dye (NQOH) 
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Fic. 2. Variation of optical density at A=600 my with pH of the solution. 


which acts presumably as an ampholytic substance exhibiting the following 
equilibria : 


H.O 
-HO + tHNQOH z NQOH * NQO- + H+ 
a 

It can be shown that at half-transformation points of the inflexions in 
the curve given in Fig. 2, the pH is roughly equivalent to pK of the system; 
this leads to the deduction from data in Fig. 2 that pKp and pK, are roughly 
7-3. and 8-1. It is to be noted that the system is complicated not only by the 
two dissociation processes but also by close proximity of the values of pK, 
and pKp, indicating that at any pH value, the optical density of the system is 
contributed not only by the cation or anion and neutral molecule, but by all 
the three forms of the substance. 


The following considerations employed for calculation of the optical 
density at any pH value based on the values of 7-3 and 8-1 for pKy and pKg 
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referring to the basic and acidic dissociation processes respectively, are in- 
structive: Let « be the optical density of the system at any pH value. As 
indicated above, this is contributed by the three components of the system, viz,, 
NQOH, NQO- and *HNQOH, the exact concentration of any one of these 
components being dependent upon the pH value of the solution. Thus in 
highly alkaline solutions (at pH > 10), it can be assumed that the entire sub- 
stance (relative concentration = 1) exists in the form of anion, while in strongly 
acid solutions (at pH <5) the cation form of the substance predominates; 
at any intermediate pH, it is apparent that the optical density is contributed 
by the appropriate portions, controlled by the dissociation constants Kg and Ky 
(vide infra) of the three forms of the substance. Let C,, C, and Cy be the 
concentrations of the cation, anion and uncharged components of the substance 
at any pH value; and £,, E, and Ep, represent the relative extinction coeffi- 
cients, constant characteristic quantities of the three forms, respectively, the 
optical density « at any pH is given by 


e= E,C, + E,Cy)+ EC; (1) 


Further, from the observed optical densities at pH = 4 and 12, where the 
entire substance (relative concentration = 1) is in the cation and anion forms 
respectively, E, and E, can be assumed to be 0-035 and 1-1. The value of 
E, can be computed from the relationships that 


= Co _ Ca 
pH = pKp log e* pKa + log C, (2) 


and from « = 0-8 at pH = 8-1 where C,=(C,; we have therefore from 
eq. (1) that 


« = 0:035 C, + 0°55C,+1-1C, 


The values of C,, C, and Cy, for any pH value can be computed from eq. (2) 
and from that C, + Cy+C,=1. The values of « calculated at different 
pH values are indicated by crosses in Fig. 2. It was instructive to note that 
despite the fact that the variation of the optical density with pH was not simple, 
the calculated values of the optical density agreed closely with the observed 
data. 


SUMMARY 


Absorption spectra of 4-(p-nitrophenyl-azo)-8-hydroxy-quinoline (NQOH) 
in the visible region was studied at different pH values in the range 4-12. 
NQOH exhibited two absorption maxima, one at A = 480 my in acid solution 
and the other at A= 600 mp in alkaline solution. The -absorption at 
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) = 600 mp was markedly variant with the pH of the solution; the data indi- 
cated two dissociation processes presumably due to the ampholytic character 
of NQOH 
-HO + +*HNQOH = NQOH = NQO- + Ht 
(H,0) 

with 7-3 and 8-1 for pK» and pK, referring in acidic and basic dissociation 
processes, respectively. The optical density values calculated from these 
constants agreed closely with the observed data. 
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1. INTRODUCTION 















THE Raman spectra of crystalline carbonates have been the subject of many 
investigations because of the simplicity of their structure. The carbonates 
are of special interest as they provide data necessary for the verification of 
theoretical ideas concerning the influence of crystal symmetry and the 
cation on the vibration spectra of the CO, ion and the lattice. Of the ortho- 
rhombic carbonates MCO, (with M = Ca, Sr, Ba, Pb) strontianite which 
occurs as a mineral has not been studied so far. It can be expected to yield 
some useful information concerning the vibration spectra of orthorhombic 
carbonates since the polarisability of the Sr ion is intermediate between that 
of Ca and Pb ions (Couture, 1944 and 1947). 


2. EXPERIMENTAL DETAILS 





Since strontianite is transparent to the ultraviolet, it was studied using 
the resonance radiation of mercury as exciter. A medium quartz spectro- 
graph and a Littrow quartz spectrograph with a dispersion of about 140 
and 40 wavenumbers respectively in the A 2537 region were used in the pre- 
sent investigation. The specimen employed was a lump of polycrystalline 
strontianite and as the individual crystals were not bigger than 1 mm. in 
size, only a thin layer on the surface was effective in giving the scattered radia- 
tion. Exposures of the order of an hour were sufficient to get satisfactory 
spectra with the smaller instrument while with El (large) spectrograph 24 
hours were necessary to get well exposed spectrograms. 














3. RESULTS AND DISCUSSION 


The recorded Raman spectrum reveals 13 Raman lines of which 5 may 
be attributed to the internal oscillations of the CO, ions in the unit cell and 
the remaining 8 due to lattice oscillations. The frequency shifts are 103 (4), 
115 (3), 150 (16), 184 (12), 218 (1), 238 (2), 244 (3), 254 (1), 703 (4), 1076 (20), 
1406 (-5), 1438 (1) and 1449 (1), the figures within brackets indicating roughly 
the relative intensities. The spectrum taken with the bigger instrument is 
reproduced in Fig. 1 along with its microphotometer record. Even under the 
96 


cS ee ee kl OU Oe. GP at 





id 


—7. Veo" 








Raman Spectrum of Strontianite (SrCOs3) 97 


high dispersion of the littrow spectrograph, the lines appear to be broad and 
not split up. 


The structure of strontianite (SrCO,) is identical with that of aragonite 
and belongs to the orthorhombic space group V;,2* (Pnma). The unit cell 
contains 4 molecules of SrCO;. A group theoretical analysis of the vibration 
spectra of orthorhombic carbonates has been carried out by Bhagavantam 
(1941) and by Couture (1944). The main results are summarised below. 
Of the 57 fundamental vibrations of the unit cell 12 internal oscillations and 
18 external oscillations (12 of the translational type and 6 of the rotatory 
type) are allowed in Raman effect. 


Internal Oscillatiéns—The CO; ion in the free state has a symmetry 
Dzn. Its spectrum consists of 4 distinct frequencies 14, vg, vs, v4, the last two 
being doubly degenerate. The respective frequency shifts are 1080, 850, 
700 and 1450cm.-! Because of the lower symmetry (C,) of the ion in 
strontianite the degeneracy would be removed. There may be further 
splitting arising from the coupling of the four ions in the unit cell. Contrary 
to the theoretical expectation of 12 Raman frequencies, the recorded spec- 
trum of strontianite indicates the existence of only 5 frequencies due to 
internal oscillations. It can, therefore, be inferred that the splitting arising 
from the multiplicity of groups in the unit cell is not observed in the Raman 
spectrum of strontianite. Due to the irregular nature of the specimen used, 
the recorded spectrum was not intense enough to exhibit the octave of the 
inactive frequency 850cm.-! Unlike in aragonite, this frequency is not 
observed as a fundamental in the Raman spectrum of strontianite. 


Lattice Oscillations—The observed lattice spectrum consists of 8 fre- 
quencies as compared with 14 in aragonite (Krishnan, 1950). The intense 
lines 150 cm.-! and 184cm.-! may be attributed to angular oscillations of 
the CO; groups about the principal axes of inertia contained in the plane 
of the ion (Kastler and Rousset, 1941). It is interesting to note that, while 
the line arising from the oscillations about the Y-axis is nearly the same in 
both strontianite and aragonite (150 and 154cm.' respectively), that 
arising from the oscillations about the X-axis is much different (184 and 
208 cm.-'). On the basis of the principle of equipartition of energy, it is 
known (Couture, Joc. cit.) that the ratio of the intensities of the two lines 
should be as the inverse ratio of the square of their frequency shifts. But 
it is experimentally found that I,/I; = -75 while vy?/vy?= -585. There is, 
therefore, as in aragonite, only a fair agreement and the discrepancy indi- 
cates the magnitude of influence of the crystalline field. 
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SUMMARY 


The Raman spectrum of strontianite (SrCO;) has been studied for the 
first time using the A 2537 resonance radiation of mercury as the exciter, 
== Fhe observed spectrum has been compared with that of aragonite which is 
isomorphous with strontianite. Out of the 13 frequency shifts recorded, 
* 5 have been attributed to the internal oscillations of the CO, ion and the 
* remaining 8 to the lattice oscillations. The splitting of the degenerate », 
vibration of the free CO, ion in the crystal into 3 frequencies and also the 
observed relative intensity of the lines 150 cm. and 184 cm.-? arising from 
the angular oscillations of the CO, ions about their principal axes of inertia, 
indicate a considerable distortion in the symmetry of the CO, ion in strontia- 

nite as in the case of aragonite. 
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FIG. 1. (a) Raman spectrum of Strontianite taken with a large quartz spectrograph [E, ] ; 
(6) Its microphotometer record. 
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